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We study the competition between a random potential and a commensurate potential on 
interacting fermionic and bosonic systems using a variety of methods. We focus on one dimen- 
sional interacting fermionic systems but higher dimensional bosonic and fermionic extensions, 
as well as classical equivalents are also discussed. Our methods which include bosonization, 
replica variational method, functional renormalization group (RG) and perturbation around 
the atomic limit, go beyond conventional perturbative expansions around the Luttinger liquid 
in one dimension. All these methods agree on the prediction in these systems of a phase, the 
Mott glass, intermediate between the Anderson (compressible, with a pseudogap in the optical 
conductivity) and the Mott (incompressible with a gap in the optical conductivity) Insulator. 
The Mott glass, which was unexpected from a perturbative renormalization group point of view 
has a pseudogap in the conductivity while remaining incompressible. Having derived the exis- 
tence of the Mott Glass phase in one dimension, we show qualitatively that its existence can 
also be expected in higher dimension. We discuss the relevance of this phase to experimental 
systems such as disordered classical elastic systems and dirty bosons. 
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I. INTRODUCTION 

In many systems, a competition between order and 
disorder has drastic consequences for the physical prop- 
erties. Such effects are paramount when the pure system 
has a gap in its excitations. This situation occurs in a 
variety of experimental systems. The most obvious one 
is a Mott insulator, where interactions lead to gap in 
the charge excitations. Low dimensional systems provide 
also many experimental situation where this competition 
occurs. On the theoretical side, examples include disor- 
dered spin-J chainsEJ, spin 1/2 ladders with-npn-magnetic 
impuritiescl, disordered Mott insulatorsuuMtLdoped spin 
1 chainsQ, and disordered ladder systemsDoijIId. On the 
experimental side, examples include doped spin-Peierls 
systemsE3E3, spin ladder systemstl. But such a phe- 
nomenon is not limited to fermionic system. Interact- 
ing bpaanic-.sjffitems can also lead to a Mott insulating 
phaseoOE^EZI, with which the disorder can compete. 
Using the standard analogy between d + 1 classical prob- 
lems and d quantum ones, it is easy to see that such a 
problem also encompasses elastic systems such .as. vpr.- 
tex lines in the presence of a columnar disorderE£tE£(E3. 
Other-miined elastic structures such as the charge density 
wavesEj of spin density waves for which the competition 
between a commensurate substrate and disorder easily 
occurs are also prime candidates. 

In all these systems the disorder tends to close the gap. 
In some case the mechanism is simple. Indeed when the 



ground state is degenerate and disorder lifts this degen- 
eracy an infinitesimal disorder causes the formation of 
domains and leads to a gap closureo due to the Imry- 
Ma effect. However, in most cases the ground state is 
not degenerate. In that case a finite amount of disorder 
is needed to induce gap closureO. In the latter case, the 
complete description of the gap closure is extremely diffi- 
cult with the usual analytic techniques such as the renor- 
malization group due to the absence of a weak coupling 
fixed point at which the gap would close. To address this 
problem one has to tackle simultaneously strong disorder 
and strong interactions. A question of interest is of course 
what is the nature of such transition. In particular how 
one can go from an extremely ordered (gapped) phase, 
to the (gapless) disordered one, which is known to have 
glassy properties, at least in high enough dimensions. 

Not surprisingly, given the complexity of the problem, 
very little is known. In one dimension an RG study, 
combining the RG for a pure commensurate [ ss r stemsE3'Eij 
and an incommensurate disordered systeno has been 
performed^. Since both the commensurate potential 
(Umklapp) and the disorder are relevant operators, no 
controlled analysis of the transition could be done. It 
was inferred from these studies that one goes from the 
Mott phase to the disordered one (Anderson) depending 
on which operator became relevant first. The idea of a 
direct Mott-Anderson transition seemed the most natu- 
ral one and was the one usually assumed in the literature. 
Solution on a special point (Luther-Emery line) also sup- 
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ported such conclusions! 

In the present paper we reexamine this problem. For 
fermionic systems it is of course difficult to tackle the 
interactions in general so we will mostly focus on the 
one dimensional case where the interactions can be han- 
dled via the bosonization technique. This allows us to 
derive a phase Hamiltonian that makes the connection 
between this problem and the disordered elastic prob- 
lems. A study of this phase Hamiltonian using better 
suited methods that capture some non perturbative ef- 
fects: (i) an atomic limit; (ii) a variational method; (iii) a 
functional renormalization group method, allows to reach 
the consistent conclusion that the transition between the 
Mott insulator and the Anderson phase (that we call An- 
derson glass to emphasize its glassy properties) is not di- 
rect. An intermediate phase, which has a gap in some 
of its excitations and yet is glassy does exist. We deter- 
mine the characteristics of this phase, that we call the 
Mott glass. Since the phase Hamiltonian does describe 
quantum crystals and bosons in arbitrary dimension but 
interacting fermions only in one dimension, we also give 
an excitonic argument directly for the fermions that indi- 
cates that this phase exists also for fermionic systems in 
dimensions greater than one. Some of these results were 
presented in a shorter form in Ref. 

The plan of the paper is as follows. In Sec. |l| we intro- 
duce fermionic models for t he disordered Mott systems. 
We then show in Sec. II B how in one dimension this 

to a p hase Hamilto- 
Sec. 



II C links this 



model reduces, using bosonization, 
nian that will be the core of our study, 
phase Hamiltonian with the other quantum crystals and 
classical disordered elastic systems for which our study is 
relevant. Sec. Ill is devoted to the analys is of t his phase 
Hamiltonian, using an atomic limit (Sec. Ill A), a varia- 



tional method (Sec. [II B) and a functional renormaliza- 
tion group study (Sec. [II C ). We show the existence of 
the Mott glass phase which is both incompressible and 
glassy. A reader only interested in the physical proper- 
ties of the MG phase can skip these relatively technical 
sections and go straight to the Sec. |^ where we examine 
in details the physical properties (correlations functions, 
transport etc.) Since the link between the phase Hamil- 
tonian and the interacting fermions only exists in d = 1 
we exa mine directly the fermions in higher dimensions 
in Sec. IV B and give an atomic limit argument showing 
that the physics of the Mott glass phase exists regardless 
of the dimension. Finally the conclusions can be found 
in Sec. 0. Some technical details are pushed in the ap- 
pendices of the paper. 



II. MODELS AND PHYSICAL OBSERVABLES 

A. Interacting fermions 

We want to study the competition between a Mott and 
an Anderson insulator in a dirty fermion system at com- 
mensurate filling. The prototype model for this problem 



is the extended Hubbard model with a random on-site 
potential at half-filling: 

H = -t ( c L c j> + h - c.) + Uj2 n iA n hl 

+V^2n l n 3 +^W i n il (1) 

{id) 1 

where (, ) denotes sum over nearest neighbors, a is the 
spin and and rij = + rij| is the total fermion num- 
ber on site i. Wi is the random potential at site i. For 
reasons that will become clear we also include a nearest 
neighbor repulsion V. A general discussion of the physics 
of this model will be given in section [V . Given the com- 
plexity of this model, let us first examine a much simpler 
situation in which explicit calculations can be performed. 



B. Interacting fermions in d = 1 and bosonized 
Hamiltonian 



In one dimension, many simplifications occur. Indeed 
in one dimension, it is possible to reexpress the Hamilto- 
nian ([!]) in terms of the collective charge and spin excita- 
tions of the system. This procedure is by now standard in 
one dimension and we refer the reader to Ref. p6| , p7|p^ , |29| 
for more details. In terms of the bosonic charge </> p and 
spin 4> a collective variables, (El) becomes: 
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Where p(x) is the continuum limit of the charge density 
and reads: 



p(x 



V2d x (j)p 



1 



7r (27ra) 
+p cos(V8c5 p — AkFx) 



(0) 



where p~o and is renormalized amplitude and a is a length 
of the order of the lattice spacing. All microscopic inter- 
actions are absorbed in the Luttinger parameters it p , K p . 
Spin rotation symmetry leads to = and K a = 1 at 
low energy. For very repulsive interactions (K p < 1/3) 
the Akp depaity fluctuations are the-imost relevant, as 
can be seei£3 from (||). In this limhEd, spin fluctuations 
suppress the 2kp part of the density fluctuations with 
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repect to the 4fcp. Let us note that such limit cannot 
be achieved within the pure Hubbard model with only 
on-site interactions. However, it can be achieved within 
an e xtended Hubbard model with interactions of finite 



rangetmcat^l. Thus one can study a Hamiltonian con- 
taining only the charge degrees of freedom: 



H 



dx 

— nu„ 
2tt p 



-- — / dxcos 



f H w 



(7) 



One can perform the rescaling: <f> — \/2<j)p and n = 
ILj/v2 which leads to the action, where we have intro- 



duced v — u p , K = K p 



— = ax ar , 

n J Jo \-2tvk 

9 
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(d r 0) 2 



v{d x( j>f 



■ cos 2<j) + Hw | 



(8) 



The Hamiltonian (0) also describes interacting one- 
dimensional spinless fermions in a commensurate peri- 
odic plus random potential. The lattice form of such a 
model is: 

H = -tJ2(4c t+ i +H. c.) +^2(Wi + g cos(2k F ia))c\ci 

i i 

+V^ n i n i+l (9) 



Where WiWj = DSij and kp is the Fermi wavevector of 
the spinless fermions system. In the continuum (||) leads 
to: 



H 



-ihvt 



dx 



-V I dxp(x) 2 - g I dx^ipL + ^l^r) 



dxW(x)p(x) 



(10) 



g measures the commensurate potential, V the interac- 
tion, and D the disorder strength. Upon bosonization, 
this Hamiltonian ( JlCj ) becomes (0). One can thus see 
that in one dimension, there is no essential difference in 
the charge sector between a band insulator (with a 2k f 
periodic potential) and a Mott insulator (which-can be 
viewed as a system in a periodic potentiaO). We 
stress that in order to establish the equivalence of the 
bosonized charge sector of the Hamiltonian (nj) and the 
bosonized representation of the Hamiltonian (M) does not 
rely on the equivalence of charge excitations of the Hub- 
bard model in the limit U/t — > oo with spinless fermions. 



1 Random potential 



The effects of disorder on a opa-dimensional non- 
interacting system are well knownoE3. Disorder, how- 
ever weak localizes all electronic states leading to an in- 
sulating behavior. However, contrarily to what happens 
in the case of the periodic potential, there is no gap at 
the Fermi level but a finite density of states. Also, the 
a. c. conductivity does not show a gap but a behavior 
of the form o~(u)) ~ u> 2 up to logarithmic corrections. In 
the presence of interactions, disorder can be treated by 
bosonization. For weak disorder one can separate in the 
random potential the Fourier components close to q ~ 
(forward scattering or random chemical potential) and 
q ~ 2k p (backward scattering) as: 



W[x) = fi(x) + £{x)e MkFX + C(x)e 



-l2kpX 



ill) 



and treat them separately. The Hamiltonian (Q) becomes 
for g = 



Hdis = 



dx 



Ti v 
'2tt 



(d x (f>) 



2 1 



K 



dxp(x) — d x (j>(x) 

7T 
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where fj,(x)fj,(x) = D j5(x — x') and £,{x)^*{x') — D\ ) 8{x — 
x') and D f = D b =W. p 

The random chemical potential [i can be absorbedcil in 
the quadratic part of the Hamiltonian ( |l2] ) by performing 
the transformation: 



(x) -» <f>(x) + 



K 

hv 



p(x)dx 



(13) 



Therefore, this term has no role in Anderson localiza- 
tion in the interacting system, in analogy with the non- 
interacting case£3. The backward scattering £ causes-Ap-, 
derson localization. Using the renormalization groupr 4 lr 6 l 
it is easy to see that £ is relevant for not too attractive in- 
teractions K < 3/2 and becomes of order one at a length 
scale 



7 - (^L_\ 
l °- a {l6WaK2 



f2„a \ l/(3-2*Q 



(14) 



identified as the localization length in the interacting sys- 
tem. Beyond this length Iq, the phase <f> becomes random 
and all correlations decay exponentially. 



2 Disorder and commensurate potential 

In the absence of disorder a commensurate potential 
leads to a gap opening for K < 2. When disorder is 
added to such a commensurate phase, its various Fourier 
components should be distinguished. Both the forward 
scattering and the backward scattering can compete with 
the commensurate potential, but as we have seen in the 
section above, they can lead to quite different types of 



Let us now describe qualitatively the effect of the var- 
ious components of the random potential. 
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ground state. The most interesting case is the compe- 
tition of the commensurate potential with the backward 
scattering. 

In order to understand the competition between the 
commensurate potential and the backward scattering, 
one can argue that the phase physically realized will be 
the one with the shortest correlation length. For the Mott 
phase the relevant length is the Mott length d, which is 
the inverse of the gap, or the size of a charge soliton. Thus 
if d < Iq one could expect the system to be a gapped Mott 
insulator, whereas for l < d the gap would be washed 
out by disorder and the system would be in the Anderson 
insulating phase. This qualitative argument can be put 
on a more formal basis by writing perturbative RG equa- 
tions for the coupling constant g of the commensurate 
potential and the disorder potential W . Both the pure 
commensurate case and the disordered incommensurate 
one lead to runaway flow where the coupling constant 
(resp. g and W) reach strong coupling (resp. at lengths 
d and Iq). A naive extrapolation consists in assuming 
that the phase that is physically realized is the one for 
which the coupling constant reaches strong coupling-first. 
Based on such an extrapolation of the RG analysis^ one 
thus expects a single transition between a commensurate 
(incompressible) phase and an Anderson (compressible) 
insulator. In order to go beyond this uncontrolled ex- 
trapolation to strong coupling of the RG results, we will 
use in this paper several non perturbative methods. 

Note that a complication arises from the fact that 
in the presence of the commensurate potential forward 
(resp. backward) random potential is generated by the 
backward (resp. forward) component as shown on Fig- 
ure [lj So in principle these two components of the dis- 
order should not be treated separately. However, for 
K > 3/2 the backward scattering is irrelevant so one can 
focus on forward scattering alone. On the other hand in 
the limit of strong repulsion, it can be shown that the 
closure of the gap that would be induced by a purely 
forward potential would occur at much stronger disorder 
than the one caused by a purely backward disorder (see 
Appendix |a|) . It is thus reasonable to expect that in the 
limit of strong repulsion, and for a backward disorder of 
the same order of magnitude as the forward disorder, the 
closure of the Mott gap is to be attributed to the back- 
ward component of disorder. This allows in the limit 
of strong repulsion to neglect the forward component of 
disorder altogether. Since in the following we will only 
consider this case, we will be justified in dropping the 
forward component. 

A detailed treatment of the forward scattering can be 
found in Appendix |A| 



C. Other quantum and classical elastic systems 

Besides interacting fermions in one dimension, the 
phase model describes many other physical systems both 
in one and in higher dimension. It is easiest to discuss 
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FIG. 1: The generation of effective backward (a) and for- 
ward (b) scattering from forward (resp. backward) scattering 
and commensurate potential, fi and £ denote respectively the 
forward and backward part of the random potential (see text) 
and g the commensurate potential. R and L denote respec- 
tively right and left going fermions with momentum close to 
+kF (rep. — kp). 



it in the Lagrangian path integral formulation. To fix 
notations let us first generalize the imaginary time quan- 
tum action (|J) (entering in the path integral in (d + 1) 
dimensions as J D(j)e~~ s / h ) to arbitrary dimension d as 
follows: 



, V p {<t>{x)) 



dr / d d x\ 



2nK 

V(4>(x,t),x) 



{d x <pf 



1 



2-kvK 



{d T <}>f 



(15) 



with V((j>,x)V (ft, x') = 5{x - x')R{<j) - 4>'). The bare 
disorder correlator and bare periodic potential can be 
chosen as R{4>) = 2( ^ )2 cos 20 and V p {cj>) = cos(20) 
respectively, although higher harmonics, preserving the 
7r periodicity, do appear under coarse graining and play 
an important role (even in d = 1). 

Before proceeding, and since the classical limit of this 
quantum action will be of importance below, let us note 
that it is obtained by letting h — > 0, K — > 0, K = K/h 
fixed with j3 fixed, the zero temperature limit [3 — > +oo 
(of interest here) being taken at the end (this can be seen 
e.g. by rescaling r = Tit' so as to keep the bound of 
integrations fixed as h—*Q). 

There are two types of systems which can be described 
by (H^) , quantum elastic systems with point disorder and 
classical equivalent systems with correlated disorder as we 
now describe. 
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1 Quantum crystals with point disorder 

Let us consider a quantum crystal in dimension d in 
a commensurate periodic potential plus a random po- 
tential. In that case, each particle can be described by 
its displacement with respect to its equilibrium position 
u(x), and the associated phonon modes. In general the 
displacement field has N components, with N = d for 
crystals of bosons or fermions, N = 2, d — 3 for a crys- 
tal of vortex lines etc.. At T = 0, the system can still 
have quantum fluctuations, leading to a quantum crystal 
(for a review see Ref. ^). Examples of quantum crys- 
tals to which our preseatp&tudy can apply are CDWtiJ, 
electron Wigner crystalHtJ, electrons at the surface of 
helium (d = 2), stripes in oxides. Other systems with 
N < d such as e.g. a vortex lattice at temperature low 
enough such that quantum fluctuations of vortices be- 
come important can also be studied. In that case a pe- 
riodic potential also exists from the underlying crystal 
or, for layered superconductors when the field is applied 
parallel to the layers. 

The Hamiltonian of such system can be written as: 



H — Hq + Hp + Hv 



(16) 



The harmonic part of the Hamiltonian of the system then 
reads: 



H 



d d x 



(17) 



where x = (x±, . . . , Xd), the Cy is an elastic matrix and 
the IT are the mqaienta. The particle density can be 
written as the sumEl: 



p (x) =J2S(x-R + u(R)) = p e lG ' {x ~ 



(18) 



over all reciprocal lattice vectors G, and u = Ui, . . . , Ud- 
This allows to write the periodic part of the Hamiltonian 
as: 



H = J d d xJ2v G e lG - u{x \ 



and the disorder part of the Hamiltonian: 



H = I d d x^W G {x)e lG - u(x) 

G 



(19) 



(20) 



To minimize technicalities we will not explicitly study the 
general case of an arbitrary lattice but a simpler N = 1 
version where one keeps only one component to the dis- 
placement field u. In addition of being already a good 
approximation in some cases (i.e keeping only the trans- 
verse displacement and its associated shear modulus in a 
2d lattice) the general case is rather similar up to alge- 
braic complications related to the tensor structure. Thus 
the quantum crystal with point disorder and commen- 
surate potential can be modeled by the quantum action 



( |T5| ) with the correspondence <p{x) = iru(x)/a, where a is 



the lattice spacing 

refer the reader tc 
N > 1. 



and — — the elastic coefficient. We 

2a 2 vK 

for more detailed descriptions for 



2 Equivalent classical systems with correlated disorder 

A d + 1 dimensional classical elastic system in presence 
of correlated disorder and periodic potential is described 
at temperature T c i by its partition sum: 
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dz I d d x(c(d x (f)) 2 + c i4 (d z 
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2r c / Jo 

+V p (<f>(x,z))+V(<f>{x,z),x)] 



(21) 



(22) 



where <fi(x, z) is a deformation field and L a thickness 
in the direction of correlation. For a system with inter- 
nal periodicity, such as a classical crystal or a classical 
CDW one has 2(j> = 27ra/a, a being the lattice spacing 
and u(x, z) the N = 1 displacement field. In that case the 
disorder V(4>, x) and periodic modulation (i.e the density 
for a crystal) have the same periodicity as given above. 
A prominent example is the flux line lattice in super- 
conductors (which has N — 2) in presence of columnar 
defects, c and C44 are then respectively proportional to 
the bulk (or shear) and tilt modulus. 

The two problems, i.e ( p2|) and (|l5|) are thus directly 
related via the correspondence: 



(23) 
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with K = K/h. The two equivalent models can thus 
be studied simultaneously. The classical limit of the d 
dimensional quantum model correspond to the zero tem- 
perature limit of the d + 1 equivalent classical model. 
Note that the boundary conditions may differ: periodic 
for quantum particles (antiperiodic for fermions) but usu- 
ally free for the classical system (unless artificially consid- 
ered on a torus). Note also that another correspondence 
could be defined with L — (5 and z = t/Tl. 



III. STUDY OF THE PHASE MODEL 

In this Section we study the phase model. Before em- 
barking on the heavy machinery of the replica variational 
method (in d — 1) and of the functional RG (in a d = 4— e 
expansion), we first show how the three phases of the 
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model can be obtained very simply in the following dou- 
ble limit (i) classical limit (ii) atomic limit. Perturbations 
around those limits can then be done and is not expected 
to yield drastic changes, as confirmed by more sophisti- 
cated methods below. 



A. phase diagram from the atomic limit 

We focus in this Section on the classical limit of the 
model ©, i.e % 0, K -> with K ~ K/h and (3 
fixed and further consider the zero temperature limit by 
taking (3 — * +oo at the end. We also perform the rescaling 
4> — > 4>/2 to simplify the equations. As will be shown in 
the following the phases identified here survive at small 
enough K > for d > 1. Indeed perturbations away from 
this limit are irrelevant in the RG sense. As is well known 
the classical version is still non trivial since there is still 
a competition between the commensurate potential and 
disorder on one hand and the clastic term which produces 
a non trivial classical configuration satisfying: 

= -^V^) + — (ssinOAx)) 
d(p u (x) AttK 7ra 

+v(x) cos(0°(x) - C(x))) = 0, (24) 

where we have define £(x)*/(tto) — iv{x)e l ^ x ' . There 
may be several solutions to this equation (apart from the 
global periodicity (p°(x) — » 4>°{x) + 2mir) but the phys- 
ically relevant ones that we consider here, are the ones 
with lowest energy (or action S[0 ]) which are selected 
as h -► 0+. 

We now consider the additional limit = — > called 
atomic limit because the model effectively becomes zero 
dimensional in that limit. In a second stage we describe 
the deviations from the atomic limit. We assume every- 
where that the disorder is bounded which turns out to be 
of some importance, the case of gaussian disorder being 
discussed later. 



and thus there is a single minimum (for a > 0) at </>o = 

An interesting change of behavior however occurs at 
W = g. For W < g the distribution of a is bounded 
away from with g — W<a<W + g and the new min- 
imum is distributed in the interval —<p max < C < 4>max 
with sax{4> max ) < W/g. For W > g, a is distributed in 
the interval < a < W+g and thus can take values arbi- 
trarily close to zero. Simultaneously, the new minimum 
position £' is now distributed in all of [0, 2tt}. 

Thus in this simple model two things happen simul- 
taneously as the disorder width W increases beyond 
W = g. First the distribution of the Hessian eigenvalue 
H"(<po) — a extends down to (while it is bounded away 
from zero for weaker disorder) and, second, the probabil- 
ity distribution of <f>Q changes abruptly at W > g (while 
it is bounded in a subinterval of [— tt/2, 7r/2]) below. As 
will become clear below, this abrupt change of behavior 
corresponds to a direct transition from the Mott Insula- 
tor to the Anderson Glass, which is in fact a multicritical 
point in the fi = phase diagram. 

It turns out that the above form for Hi(<j)) does not 
yield the generic behavior for d = 0. This can easily 
be seen by adding higher harmonics and we will illus- 
trate it by simply adding a small second harmonic to the 
disorder. It must be stressed that these higher harmon- 
ics are always generated in per turba tion theory beyond 
the atomic limit (see e.g. Sec. 1IIC) and that they are 
generically present in realistic models and should thus be 
included. Thus we now study: 

H 2 (cf>) = -vcos(<p-()-v 2 cos(2(f>-2( 2 ) 

- 9 cos(0) (29) 

with £2 is another random phase uniform in [0, 2n] inde- 
pendent of £. One can rewrite: 



H 2 (4>) = a(cos ip + /3 cos(2?/> - 

V = <f> - £ 

ft = v 2 /a x = C2 - C' 



2X)) 



(30) 
(31) 
(32) 



1 Atomic model: K 







Dropping the elastic term in (g^) we are thus left to 
study the d — model Hamiltonian in the classical limit 



Hi{(f>) = — ucos(0 — C) — gcos((f>) 



(25) 



with C a random phase distributed uniformly in [0, 27r] 
and v a random variable, which we choose to have a 
bounded support — W < v < W. In the absence of disor- 
der the minima are at <fi — 2ixn. In the presence of disor- 
der this model has (up to the periodicity 4> — > 4> + 2nn) a 
unique local (and global) minimum with probability one 
for all parameters. Indeed one can rewrite: 



Hx{(j)) = -acos(4>-C) 
a = \Jv 2 + 2ugcos(£) + g 2 



g + ve 



-iC 



(26) 
(27) 
(28) 



and a, (' as above. Note that \ is still uniform in [0, 2ir]. 
Let us consider a fixed f3. It is easy to see that for (3 < 1/4 
there is a unique minimum for any value of \ an d the 
situation is similar to the one discussed above. But, as 
soon as (3 > 1/4 there is a value of x f° r which there is 
a second minimum. At (3 — l/4 + the second minimum 
appears for ^ = 0. Thus, new things happen in this 
model. The phase diagram is shown in Fig. ||. Let us 
consider for simplicity the case where v\ = v and v 2 are 
fixed and positive and only the phase are random (the 
general case is similar). 

There are three phases in this simple, exactly solvable, 
model: 

(i) For v\ < 5, a is bounded from below (a > g — V\) 
and for small enough v 2 < \{g — v\) a single minimum 
exists: this corresponds to the MI phase as shown in 
Fig. I and Fig. |. 

(ii) For v\ < g and v 2 > j(g — v±) two minima exist. 
This corresponds to the MG phase. Just above the line 
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FIG. 2: Phase diagram in the atomic limit: V\ and «2 are 
the disorder strengths of the two disorder harmonics and g 
the strength of the commensurate potential. The three dif- 
ferent phases are: MI which has a unique local minima with 
probability 1; MG which corresponds to a non zero probabil- 
ity to have two local minima and zero probability that any 
minima lie outside of ] — tt/2, tt/2[; AG which has a non zero 
probability to have two local minima and finite probability 
that minima are outside of ] — tt/2, tt/2\ (and thus that there 
are kinks, see text). 

v 2 = \{g — v\) the equilibrium position is ip — x ~ 0, 
4>o ~ C and |0 O | < 4> m < n/2 with sin0 m = v\/g. Thus, 
with probability exactly one the two minima remain in 
the wells of the original cosine, i.e the probability that 
the new minimum is outside of the interval [— 7r/2,7r/2] 
is exactly zero. This is represented in Fig. ||. 

(hi) For Vi > g, a has a finite probability to be arbi- 
trarily close to zero and one easily sees that the prob- 
ability of having two minima is nonzero and the prob- 
ability that the new minimum is outside of the interval 
[— 7r/2, 7r/2] is non zero. This phase corresponds to the 
Anderson glass as is shown on Fig. ^. 

2 Expansion around the atomic limit, d > 1 

We can now expand around the atomic limit and con- 
sider large but finite K. Let us consider d = 1 for simplic- 
ity, but similar arguments apply to any d > 1. We must 
construct the classical configuration of the Hamiltonian: 

H(4>) = J !(V s 2 -«coB(0(a;)-C(aO) 

— v% cos(2(f)(x) — 2( 2 (xj) — gcos((j)(x)) (33) 

with c = v/ttK. Assume again for simplicity v and v 2 
constants. Let us work in the limit of elastic coefficient 
c very small (large but finite K or equivalcntly of a long 
correlation length for the independent random phases 
£(x) and ( 2 (x)). Then we can think of the model as a 
succession along x of d = models (slices) with different 
realization of the disorder, and consider e.g. a discretized 




FIG. 3: Pure Mott Insulator phase without disorder. A 
ground state in the classical limit (f>o(x) = is represented 

version: 

H (<t>) = ^^n+l-^n) 2 ~VCOs((t) n -Cn) 
n 

-v 2 cos(2</>„ - 2C 2 ,„) - g cos(<^) (34) 

where the and £2,™ are independent from slice to slice. 

Let us think of a formal perturbation in the elastic coef- 
ficient c. For c = we know the minima for each slice, an- 
alyzed above, noted 0q + 2uk n . For c > 0, one easily sees 
that to lowest (naive) order in c, to construct the minimal 
energy configuration one must first choose which minima 
of successive slices to connect together. Small shift in 
minima positions and the ensuing changes in minima en- 
ergies is formally of higher order (they become relevant in 
the MG and AG phase and are discussed below but they 
do not change our main argument here). Thus, at each 
x slice we must choose whether to connect </>q with <j)^ +1 
(defined as the minimum in [— it, it] of the corresponding 
d = model) or to 0q +1 ±27t (which we call a kink). This 
is equivalent (formally to lowest order in c) to choosing 
iteratively the one of the three minima which minimize 
the distance min m= o,i,-i(|0o — 0o + 2mir\). The way 
these minima are connected define the three phases. 

(i) In the MI regime defined above («i < g, v 2 < (g — 
Vi)/4) connecting the minima is obvious and leads to 
a typical ground state as represented in Fig. ^, which 
defines the MI phase. 

(ii) In the MG regime (vi < g, v 2 > (g — «i)/4) 
the key observation is that connecting the minima re- 
mains unambiguous thanks to the fact that the distribu- 
tion of the position of these minima remains confined 
with probability 1 within [—<j) m , <j) m ] with <j) m < tt/2. 
Indeed as long as 4> m < tt/2 the distance for a kink 
\<J)q ~ 4>q +1 ±2tt\ > 2tt — 2<ft m is always larger than 2<p m 
the distance to connect two minima maximally separated. 
Thus the elastic energy always penalizes kinks. One thus 
finds that the string 0q will remain confined, for small 
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FIG. 4: Mott Insulator phase with disorder. The ground 
state in the classical limit (f>o(x) is represented. It is only 
slightly deformed with respect to 4>o(x) — 0. No other local 
minima exists (up to the global periodicity 4>o(x) — > (f>o{x) + 
2tt). The stability eigenvalues of the Hessian matrix at (f>o(x) 
are strictly positive and bounded from below by a positive 
number m 2 R . There is a gap ~ m\ in the conductivity. The 
compressibility is still zero since the response to a tilt hd x (f) 
vanishes as no kinks exists 



c, within the interval [— <f/ m , (f)' m ] with <j)' m < n/2. This is 
represented in Fig. || and corresponds to the Mott Glass 
phase. 

(hi) In the AG regime (vi < g), in each d = slice 
4>q can be in any position of [— 7r, 7r]. Thus connecting 
the minima becomes qualitatively (although not quanti- 
tatively) the same problem as for the standard CDW 
without the periodic potential. This is the Anderson 
Glass phase. 

Thus these semi-rigorous arguments show that the 
three phases which already exist in the atomic limit do 
survive upon adding a small c > 0. This also implies 
that for very small c we expect the phase boundaries to 
be close to the ones at c — (i.e of the d — atomic phase 
diagram Fig. ^. In practice there are of course some limi- 
tations. The perturbation in c can be used, strictly, only 
in the MI. In the MG perturbation in c cannot strictly 
apply as the effect of c on 4>q is of order c/a and thus 
generate terms as c 2 /a in the energy, which dominate 
over c for small a. Thus, as soon as the Hessian eigen- 
values extend to perturbation theory fails. In the MG 
this just means that because of possible multiple minima 
the ground state cannot be determined perturbatively. 
This however only concerns the precise position of </>° in 
the well, but does not have any consequence on the fact 
that all </>° remain in a single well. These effects which 
go beyond perturbation theory will change the precise 
dependence of the boundaries as a function of c from 
a naive perturbative estimate, but we believe that the 
boundaries are continuous as c — ► 0. The d = phase 
diagram should thus also give approximately the d = 1 



5(j) finite 




FIG. 5: Mott Glass phase: the ground state in the classical 
limit 4>q (x) is represented. Other metastable states exists, and 
the stability Hessian matrix at 4>o(x) spectrum extends down 
to 0. There is no gap in the conductivity. The wandering 
of the ground state along x is bounded as \4>o(x) — <f>o(x')\ 
is finite (and the x averaged positions is at <j> — 0). The 
compressibility is still zero since the response to a tilt hd x (j> 
(i.e a change in chemical potential) vanishes as no kinks exists 
between the well separated original minima of the cosine (j> = 
2nir. 



one for small c. 

A useful approximation, used in Sec. Ill C| to study the 
MG/AG transition and which consists in replacing the in- 
teraction term by a mass ^m 2 (j) 2 , can be checked in the 
atomic limit. It gives correctly the multicritical point at 



v = m since it has a transition for W — m (for W < m 
there is a unique minimum with probability 1 while for 
W > m 2 there is a finite probability that there is a sec- 
ond minimum) . By construction this approximate model 
cannot distinguish between the MG and the AG phase, 
and is only used to describe the MI/MG transition. 

One must emphasize that none of the transition in the 
d = model in Fig. || survives if the local distribution of 
disorder is gaussian. Indeed the probability of a second 
minimum is always non zero, although it can be exponen- 
tially small, resulting in sharp crossover behavior rather 
than transitions. It is however not so clear whether the 
same apply in higher dimension. We know that the gap 
is robust to small bounded disorder. However for un- 
bounded disorder, formation of terraces if size L become 
energetically favorable in rare regions with exponentially 
small probability. This leads to a gapless spectrum, with 
an exponentially (in d > 1) small density of state at low 
energy. Note that although the gap itself cannot be used 
anymore as an order parameter one can still clearly de- 
fine a phase transition between MI and MG since in the 
MG phase the spectrum is expected to become algebraic. 

To conclude this Section we have established that for 
bounded disorder three phases (MI, MG, AG) already 
exist in the atomic classical limit. Previous attempts 
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FIG. 6: Anderson Glass phase: the ground state in the clas- 
sical limit 4>o(x) is represented. Many other metastable states 
exist, and the stability Hessian matrix spectrum extends down 
to 0. There is no gap in the conductivity. The wandering of 
the ground state along x, \4>o(x) — 4>o{x')\ is unbounded. The 
compressibility is non zero as the ground state reorganize in 
response to a tilt hd x (j> 0- e a change in chemical potential) 
because kinks of energies arbitrary close to zero now exist. 



at analyzing the classical limited assumed that beyond a 
length Lq the distribution of the phase becomes random. 
As we find here, this is incorrect for weak disorder and the 
phase has instead a narrow distribution around (f> = 0. 
Thus these arguments missed the existence of the Mott 
Glass phase and the correlation length Lq identified in 
Ref. fry is not the correct one. In order to get more 
detailed information on the three phases in d = 1 and 
higher we now turn to more sophisticated methods. 



B. Replica Variational Method 

Let us now use a Gaussian variational method@@ to 
study the action originating from the Hamiltonian (|l 2[) . 



We only retain in this section the backward scattering 
which is the one leading to localization. The effect of 
forward scattering is analyzed in appendix |A]. Let us 
recall that backward scattering is generated even if only 
the commensurate and forward scattering are present. 



1 Derivation of the variational equations 

Let us first recall briefly the principle of method we use. 
The partition function of a disordered quantum system 
is 



Z = J d[<f>]e- 



.Slo] 



(35) 



where S[4>] is the Euclidean action that depends explic- 
itly on the quenched disorder. In physical problems, one 
usually needs the average of free energy 



F 



1- 

y 



-InZ 



(36) 



where the overbar denotes a disorder .average. This av- 
erage can be done via the replica tricke3: 



In Z = lim 

n— >0 n 



One has for integer n: 



n <w 



n 



(37) 



(38) 



where S rcp . is a disorder free quantity that depends on 
the n fields <fi a - In the end, one has traded the disor- 
der average to a analytical continuation from an integer 
number of fields n to n = 0. In our particular case, after 
having averaged over disorder the action (g|) leads to the 
replicated action: 



J 



E 



dxdr 
2nK 



W 



(2irah) 2 



v{d x (f) a f 



(d T 4>a) 



irah 



dxdr cos 2(p a 



n ,/3 Pf3 

V] / dx I dr dr' cos(2((f) a (x,T) - cj> h (x,T'))) 
,J Jo Jo 



(39) 



in which one has to take the limit n — > 0. One way 
to perform this limit is to use a Gaussian Variational 
Method (GVM) initially introduced to study classical dis- 
ordered systems-|Such as random het erap olymersLj, ran- 
dom manifoldsc3, and vortex Iatticesc3c3. This Ansatz 



has been extended to treat correlated disorder and thus 
to apply to quantum systems as wellEHl. Since this method 
has been shown to describe with a good accuracy both 
the pure commensurate phase and the Anderson insula- 
tor, on can expect to get also good results in this more 
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complicated situation. 

The ansatz consists in finding the "best" quadratic ac- 
tion 

So= 2f3fiYl J 2^M<l,u n )G-^(q,uj n )(j) b (q,uj n ) (40) 



associated with 



with : 

n-if \ (M 2 +^ 2 ) x f \ 

vG ab il,^) = o ab - a ab (q,u)) 

i.e. the one that minimizes the trial free energy: 



-Fvar = F + — (S Iep . — 5*0)50 



(41) 



(42) 



Where (. . .)s designates the averages performed with 
respect to the Gaussian action and Fq is the free energy 



(43) 



In this method the full Green's functions G(q,ui) are 
the variational parameters. The derivation of the saddle 
point equations is performed in appendix [b]. Introduc- 
ing G c (q,uj n ) = J2b G ab(q,u n ), B ab (x,r) = G aa {x,T) - 
G a b{x, t) and parameterizing these nxn hierarchical ma- 
trices using Parisi's ansatz (G a ^b — > G(u), G aa — > G, and 
similarly for any quantit y s ee appendix g), we obtain the 
saddle point equations (|B6|): 



J 



G~ 1 (q,u> n ) = A(V + — ) + — exp (-27iG(z = 0, 



h 
2W 



0) 



dr(l — cos(cj„t)) 



exp(-4KB(x = 0,r)) 
I 



du exp(—ATiB(u)) 



(44) 



with 



a(q, u> n 



2Wv 
(ira) 2 



/3exp{-MB{u))6 Un , (45) 



In the absence of the commensurate potential, ( f44| ) are 
known to lead to a one step replica symmetry .broken 
solution describing the Anderson localized phasec3, the 
replica symmetric solution being unstable. On the other 
hand, in the absence of disorder the commensurate phase 
(i. e. the Mott Insulator) is obviously replica symmetric. 
Therefore, we will search first for a replica symmetric 
solution that we expec t to be associated with a Mott 
insulating phase in Sec. [II B 2. Above a certain disorder 



this solution will become unstable, and one has to turn 
to replica symmetry broken solutions. As we will see in 
Sec. IIIB3, in the presence of the commensurate term, 



there is, besides the saddle point solution corresponding 
to the Anderson insulator, there is room for a third saddle 
point solution corresponding to the Mott glass phase. 



2 Gapped replica symmetric solution: Mott Insulator 

For the replica symmetric solution, G(q,uj ni u) 
G(q,u> n ). The saddle point equations then read: 



vG c 1 (q,u n ) 
G(q,uj n ,u) 



— ((vqf+ul)+m 2 +I{u) n ) 

2W/3K 2 e- ttG °""i M „ |0 

(ha) 2 ((vq) 2 + irKm 2 ) 2 
jg^ e -2ftG(0,0) 

2^L e -«G.(o.o) f h dT(e« G -<«W - 1)(1 - cosKr)) 



(46) 
(47) 
(48) 
(49) 



r 



In the general case, the equations (^)-(^) can only be 
solved numerically. However, in the limit h , K ^ 0, 
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K = K/h fixed, it is possible to solve these equations 
analytically. This is due to the fact that in that limit to 2 
has no dependence on I(u) n ), so that to is given by the 
simple equation: 



Agv 



■ exp 



WK 1/2 



y2 7r 3/2yl/2 m 3 



(50) 



The self-consistent equation ( p0| ) has always the trivial 
solution m — 0. Let us determine under which conditions 
it can also have a non-trivial solution with ra/0. In or- 
der to obtain the answer in physical terms it is convenient 
to reexpress all quantities as a function of the physical 
lengths Iq and d: 



1 16WK 2 



'n 



(aw) 



1 

7r- 



*gK 

(av) 



(51) 
(52) 



They correspond respectively to the localization (pin- 
ning) length in the absence of commensurability, and to 
the soliton size of the pure gap phase. We introduce the 
length 



(■nKm 2 ) 



(53) 



which, as we will see, is the correlation length in the pres- 
ence of both the commensurate potential and the disor- 
der. One can then rewrite ([sb"|) as: 



1 



1 

d 2 



exp 



16 Wo 



(54) 



1 /3 

For Zo/rf>^(x) ' this equation admits two solutions. 
It can be seen by considering the limit of large disorder 

1/3 

that solutions with Z /£ < \ (I) are spurious. Thus, 
for l /d > i f^J , we have a unique solution of (54). 

1/3 

For lo/d < i (-2^) , there is no solution, which means 
that the replica symmetric solution with a mass is unsta- 
ble. Since it is knownE2l that a replica symmetric phase 
with no mass is unphysical for h — > 0, this means that 
for large enough disorder, we obtain a breaking of replica 
symmetry. 

Let us now examine the equation for I(u> n ). An ex- 
pansion around h — in equation ( ft9| ) gives the self- 
consistent equation for I(uj n ): 



8Wv 



G c (x = 0, t)(1 — cos(uj n T))dT (55) 



Going to Fourier space and performing the q integration 
leads to to the final form: 



AWvK 



1 



1 



w2+rf(m 2 + %)) 



(56) 



/(CO) 

h=2- 



CO/CO* 



FIG. 7: The function / as a function of ui for A = 1 (full 
curve) and A = 2 (dashed curve). Note that as A = 2, / 
starts linearly with frequency, contrarily to A < 2 for which / 
is quadratic at small u. 



Obviously, 

\ V nKm J 
where / satisfies the equation: 



f(x) = A 



fix) 



and 



4WK 1/2 v 



r 3/2 Q ,2 m 3 4 Wq 



1 



(57) 



(58) 



(59) 



As can be seen from ( flq ) , m defines the correlation length 
£ in the presence of both the commensurate and the ran- 
dom potential. 

Once to 7^ is known from (|50|), A and therefore / 
and I(uj n ) are entirely determined. The above equations 
thus completely fix all the parameters of the gapped RS 
phase. 

For A < 2, there is a physical solution of ( |58l) such 
that \im x ^± oc f(x) = 1 + A and for x <C 1 , f(x) 
1 + ax 2 + o(x 2 ) with a = A/ (2 — A). The corresponding 
behavior of f(x) as a function of u> is shown on Fig. [?] for 
A = 1. For A = 2, for x < 1, f(x) = l + ^=\x \ +0{x 2 ) 
and lim^^oo f(x) = 3. The corresponding graph of f(x) 
is also shown on Figure ^. For A > 2 ( |58| ) has no physical 
solution. Thus, A = 2 is the boundary for the gapped 
RS phase. Putting A = 2 in (pj) leads to l /£, = \ > 

\ (f) 1 ^ 3 - Reinjecting this value in (E3) gives 



lo 
d 



1 i 



(60) 



This point is in the domain where (54) still has solutions 



Therefore, as disorder increases the system attains the 
point where A = 2 &e/ore reaching the point where m = 0. 
Beyond the point A — 2 the replica symmetric solution 
becomes unstable. This leads us to consider a replica 
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symmetry breaking solution of Eq. (Q) for ^ < e^/2 
allowing for a non-zero m. The corresponding phase will 
thus not be the simple Anderson localized phase expected 
from the strong coupling RG argument 62 . 



3 Replica symmetry broken solution 

In the preceding section, we have seen that in the limit 
fi->0a replica symmetric solution of equations (|44|) can 
exist for A < 2 but is unstable for A > 2. In this section, 
we consider in the limit h — > the one-step replica sym- 
metry breaking (RSB) solution of equations (44j). Such 
RSB solution should be valid for A > 2, and is known to 
correctly describe the simple Anderson insulating phase 
for the simple disordered caseEj. Compared to the case 
in the absence of commensurate potential we have to al- 
low, in our RSB solution, for m 2 ^ 0. However, con- 
trarily to the RS—case, a RSB solution with m 2 = is 
perfectly possibleo and corresponds to the case in which 
the commensurate potential is completely washed out by 



the random potential. 

Two scenarios could thus be a priori possible. Either as 
soon as the replica symmetric solution becomes unstable 
one obtains a RSB solution with m 2 = similar to the so- 
lution of Ref . |2^ or there can exist an intermediate regime 
with both RSB and m 2 ^ 0. The first case would corre- 
spond to the simple scenario, suggested from extrapolat- 
ing the RG of a direct transition between the commen- 
surate phase and the Anderson insulator. On the other 
hand the behavior in the RS solution strongly suggests 
the existence of an intermediate phase, the Mott Glass. 
Indeed in the RS phase the optical gap closes, leading to 
a conductivity much similar to the simple Anderson one, 
whereas the compressibility remains zero. This suggests 
that all effects of the commensurate potential have not 
yet disappeared, and that the system is not in a sim- 
ple Anderson regime. As we will see, this is this second 
possibility that obtains, leading thus to a much richer 
physical behavior than could have been guessed from the 
RG extrapolations. 

The saddle point equations (jij) are first rewritten as: 



vG-^q^n) = -={{vqf +ujl)+m 2 + Ei(l - <J„. ) + J(w„) 
irK 



Si 
a{u) 



2Wv 



2 +£ .2. 



e -ihB(r) _ e -4hB(u>u c ) 



{-Kaffi J 

u c {a(u > u c ) — a(u < u c )) = [a]{u > u c ) 
2Wv 



(1 — COs(iU n T))dT 



(7ra) s 



-hiB(u) 



-4hG{0) 



(61) 

(62) 
(63) 
(64) 

(65) 



where we look for a one step RSB solution, as is adapted 
to d = 1+1. The parameters m, Si and the breakpoint u c 
have to be determined self-consistently. The full solution 
is given in Appendix |^. The parameter u c is determined 
from the maxginaiity of the repiicon condition, which has 
been shownEj to be the correct condition to impose. This 
leads to I(u> n ) oc | Lu n \. One can check (see Appendix ^ 



and Sec. [II B 2 ) that this condition is also satisfied by 
I(io n ) at the limit of stability of the RS solution. The 
two other parameters m and Si depend on the ratio d/lo. 

As is shown in Appendix ^ for d/lo > 1 .86 one has m = 
and Si ^ 0. This solution is thus similar to the one of 
a disordered system without the commensurate potential 
and corresponds to an Anderson glass (insulatacj). Such 
a phase has no gap in the optical conductivity. Since 
there is no mass in the propagator ( |6l| ) for u> n = the 
compressibility is finite. 

On the other hand for d/lo < 1.86 the solution has a 
finite mass. This regime is thus different from the sim- 
ple Anderson insulator. The physical properties (con- 



ductivity, phason density of states, compressibility) will 
be discussed at length in Sec. |y|. An important result 
of Sec. IV is that because of the presence of the mass 



the system is still incompressible while having the con- 
ductivity ofLan Anderson Insulator. This is the Mott 
glass phasec3 which shares some properties of the Mott 
Insulator (incompressibility) with those of a glassy phase 
(breaking of replica symmetry). 

Therefore, the physical picture is the following: for 
d/lo < 2e -1 / 4 , one has a replica symmetric solution 
with a gap in the conductivity, the Mott Insulator, for 
2e -1 / 4 < d/l < 1.86, one has a RSB solution without a 
gap in the conductivity but zero compressibility, the Mott 
Glass and finally, for d/lo > 1.86, there is a finite com- 
pressibility and no gap in the conductivity, the Anderson 
Insulator. In other words one recovers the solution of 
Ref. ^ not as soon as d/l > 2e~ 1 / 4 as we would expect 
from extrapolations of the perturbative d = 1 renormal- 
ization group calculations, but only at the higher value 
d/lo > 1.86. This is due to the formation of an interme- 



13 



Replica 


Replica Symmetry 


Replica Symmetry 


Symmetric 


Breaking 


Breaking 


Gap in 
conductivity 


2 

a(CO)ocQ) 


a(Q))°cCQ 2 


Incompressible 


Incompressible 


Compressible 
1 ► 



U.558 



11.86 



d/lo 



FIG. 8: The phase diagram of the commensurate system 
with backward scattering only as a function of d/lo (disorder 
increases when d/lo increases). At weak disorder, d/lo *C 1 one 
obtains an incompressible phase with a gap in the conductiv- 
ity i.e. a Mott Insulator. For strong disorder, d/lo 3> 1, one 
has a compressible phase with a conductivity that behaves as 
(j{uj) oc u 2 i.e. an Anderson glass. The surprising feature of 
this phase diagram is the appearance of an intermediate in- 
compressible phase (like the Mott insulator) having the same 
conductivity as an Anderson glass for d/lo ~ 1- 



diate phase, which is both incompressible but without a 
gap in the conductivity. This intermediate phase being 
an intermediate coupling one, the failure of the pertur- 
bative renormalization group to predict its existence is 
not a surprise. In a forthcoming section, we will discuss 
its properties in detail. The phase diagram as a function 
of d/lo is represented on Fig. |8[ 

Let us remark that all transitions appear to be first 
order within the GVM formalism. 



C. Mott Insulator to Mott Glass transition: 
Functional Renormalization Group study and 
classical equivalent 

In this Section we study the phase model (|l5|) in arbi- 
trary dimension d using a renormalization group method 
perturbatively controlled in d — 4 — e and small Ti. It 
provides useful information on interacting fermions with 
disorder by continuation down to d = 1 (as we do not 
expect drastic changes in this model down to d = 1 for 
small K) . We will perform the analysis in the notations of 
the classical equivalent model, but also give some conclu- 
sions in terms of the parameters of the quantum model, 
via the relations (p3f). 

As is well known for class i ca l problems such as man- 
ifolds in random mediac3c3l£], the Functional Renor- 
malization Group (FRG) provides an alternative to the 
Gaussian Variational Method. The FRG treats accu- 
rately the non linearities and does not use replica sym- 
metry breaking. When the two methods are supposed to 
be exact and can be compared they, do agree, as found 
for the random manifold probleaao for N — > oo, and 
generally give consistent physicsc£l. Here we will use the 
FRG as a check of the correctness of the prediction by the 
GVM of the Mott Glass phase, and as a way to obtain 
additional detailed information on the Mott Insulator to 



Mott Glass transition. 

Although this is not the route we follow here one can 
apply the FRG method directly to the model @-(|2|). 
This am ount s-, to generalize to correlated disorder the 
study ofiaaa made for the case of uncorrclated dis- 
order. It shows that a commensurate potential becomes 
relevant and may lead to a description of the transition 
between the gapless Anderson Glass to a gapped phase. 
But this route would not give us the information we want 
about the nature of the gapped phase, i.e it cannot dis- 
tinguish between the M©itt-iplass and a Mott Insulator. 
Indeed, the approach oe3Gj£J fails to describe the phase 
where the commensurate potential is relevant since the 
FRG then flows to strong coupling. 

Thus, in order to test for the existence of the Mott 
Glass we consider from the start a situation where the 
commensurate potential is relevant and replace the full 
model (p~5|) — ( 22 ) by an effective model in which the sine 
Gordon term is replaced by a quadratic mass term: 



(66) 



This should be a reasonable approximation when the sine 
Gordon term is relevant, and is in the spirit of the Self 
Consistent Harmonic Approximation (SCHA). Our ap- 
proximation amounts to neglecting some soliton excita- 
tions by giving them large energy and to neglecting the 
renormalization of the gap by disorder. This simplified 
model has the merit to be amenable to a perturbatively 
controlled study in d — 4 — e. We will show that it does 
exhibit a phase transition at T c \ = (K = 0) which sur- 
vives at T c i > (h > 0) and can be identified with the 
Mott Insulator to Mott Glass transition. This model thus 
allow to study the formation of the Mott Glass. 

We have obtained the FRG equations for the effective 
model (|6^). They can be derived by integijated|-Qi 
eratively short wavelength modes, extendingc3c3'E3'l 
They are obtained in terms of the running dimensionless 
disorder A(0, 1) = —R"(<f>,l), the running dimensionless 
temperature Tj ~ T, both defined in the Appendix 
and the tilt modulus 044(1). These RG equations read: 

d,A(0) = eA(^)+f,A"(0) 

-/(O(A"(0)(A(O) - A(cj>)) - A'(0) 2 X67) 
diC44 = -/(0A"(0)c 44 (68) 



with e 
and 



4 — d, setting c = 1 (as c is not renormalized) 
1 



f(l) 



(1 + fie 21 ) 2 



(69) 



comes from the integration of the high momentum 
modes. Here \i = (ma) 2 . These FRG equations are an- 
alyzed in Appendix [5]. We describe here only the main 
results. 

At T — (i.e h — ► for the quantum problem) we 
find that there is a phase transition. One can measure 
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FIG. 9: Phase diagram of the effective model from the FRG, 
at T = (ft = 0). J? c ~ (l/W) 1/(4_d) is the Larkin length (lo- 
calization length) and parameterizes disorder strength. These 
phases and transition survive at T > (K > 0). 



the strength of the bare disorder using the Larkin length 
R c ~ (l/Af =0 («/> = 0)) 1 /( 4 -<0 of the problem without 
commensurate potential (corresponding to the localiza- 
tion length for the d = 1 fcrmion problem). Then we 
find that for a given m, there is a transition at a critical 
disorder strength 



Rc 



1 



l/(4-rf) 



AL o «> = 0) 



c- 



1 



(70) 



where C is a constant. The two phases are: (i) For strong 
disorder R c < C± we find that A 2 (Z) = -A"(0,l), the 
fourth derivative of the renormalized disorder correlator 
Ri(u), becomes infinite at a finite scale R*(m), i.e the 
disorder correlator becomes non analytic and develops a 
cusp singularity at a scale R*(m). For the problem in 
the absence of a mass the cusp generation at the Larkin 
length R* c {m = 0) = R c is well known to be associated 
to the existence of many metastable states beyond R c . 
This cusp generation is associated with the apparition of 
the transverse Meissner effect in vortex lattices pinned 
by columnar disorderE2L(and to the appearance of RSB 
in the GVM treatmentCj). As will be discussed below 
this phase corresponds to the Mott Glass; (ii) For weak 
disorder R c > the flow is cut by the presence of 
the mass before a cusp can be generated. This phase 
does not exhibit metastable states and corresponds to 
the Mott Insulator. 

Since the mass can be chosen arbitrarily small the 
study is thus perturbative in disorder in d = 4 — e for 
the model (f36|). It is interesting to note that this T = 
transition exists both for correlated and uncorrelated dis- 
order. However, this transition is stable to finite temper- 
ature only for correlated disorder. Indee d for point like 
disorder the temperature rounds the cuspr 5 lP 6 l which im- 
plies that there can exist no sharp distinction between 



the two phases at finite temperature. In addition the 
quadratic part of the Hamiltonian is not renormalized 
by disorder, and thus even at T = there cannot be 
any signature of the transition on two point correlation 
functions of <j>. Thus it is possible that the transition ob- 
served in Ref. (s?] is an artefact of the method used. On 
the contrary for correlated disorder, there is a genuine 
transition and because of the lack of rotational invari- 
ance (in (x, r)), the existence of the cusp and the tran- 
sition directly affects for correlated disorder, two point 
correlation functions. 

The FRG gives immediate information on the renor- 
malized tilt modulus (see ( |23| ) C44 = 044(4-00). Since this 
is also the coefficient of oj 2 in the Green function (c/xp) one 
can infer from that that if C44 is finite the Green function 
is likely to remain analytic and thus that there is a gap in 
the conductivity. If C44 becomes infinite then the Green 
function is not analytic and no gap should exist in the 
conductivity. The FRG gives that: 



c 44 (+oo) 



'Rc 



-) 



1 



C44 (0) (Rcy f EM J 



(71) 



and thus we find that the phase (i) above which corre- 
sponds to a cusp (C44 = 4-00 ) can be identified with the 
Mott Glass, while the phase (ii) above which corresponds 
to no cusp ( cf 4 < +00) can be identified with the Mott 
Insulator. The gap itself can be estimated as: 



A = to/ 



(Rc--) 

TO 



(72) 



thus it vanishes linearly A ~ (R c - at the MI to MG 
transition. Correlation functions are also estimated in 
the Appendix 

One of the most crucial test is to show that the above 
MI-MG transition survives to quantum fluctuations h > 
(thermal fluctuations T > for the d + 1 classical 
model) . This is the case and the calculation is detailed in 
Appendix There is no doubt that the MI survives to 
quantum fluctuations, however it was less obvious that 
the MG would survive. Indeed the cusp i s round ed by the 
effective temperature variable TJ ~ T/yJ cn{l). However 
the key point is that ca(l) becomes very large as soon as 
the second derivative A2 grows, and as a consequence the 
effective temperature T renormalizes to exactly zero at a 
finite scale (as it does in the absence of a mass). Thus 
the Mott Glass survives at finite temperature. A similar 
phenomenon was also found recently in the dynamics of 
classical periodic systems with correlated disorderEl 

To conclude this Section, the FRG shows, within a 
d = 4 — e analysis of the effective model with mass, that a 
transition exists at large K in the quantum problem (and 
at low temperature in the equivalent classical problem) 
between a MI phase at weak disorder with analytic Green 
function, no metastable states and a gap in the conduc- 
tivity and a Mott Glass phase with metastable states, no 
gap in the conductivity at stronger disorder. It allows to 
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predict that the conductivity gap should close linearly at 
the transition (at least in the limit of small K — > 0) . 



IV. PHYSICAL PROPERTIES, RESULTS IN 
D = 1 AND EXTENSIONS TO HIGHER 
DIMENSIONS 



the phason density of states of the system. One of the 
most striking difference between the Anderson Insulator 
(AI) and the Mott Insulator (MI) is that the former is 
compressible whereas the latter is incompressible. The 
compressibility is given, in any dimension, by 



A. d= 1 

1 Compressibility , density of states and correlation 
functions 

In this section, we define and calculate equilibrium 
thermodynamic quantities such as the compressibility or 



J 



i r rP n 

C (q,u n ) = - d d x dre^ x -^ T \T T (n(x,T) - (n(x, r)))(n(0, 0) - (n(0,0)»> 



(73) 



where n is the density. This leads to the average static 
compressibility Xs — hm 9 ^o(li m w ^o X{Qi UJ ))- In d = 1, 
using the bosonic expression for the density ( f73| ) leads to 



X S = lim lima G c {q,u>) 

q— >0 uj—fO 



where 



J q,ioH J —q, — uj / 



(74) 



(75) 



Another thermodynamic quantity of interest is the 
phason density of states: 



p(u) = $s[iu} n G(x,X,iLO n )] |ia>„->w+i0 4 



7TC 



(76) 



In the K — ► limit, ( |76[ ) describes the phason density of 
states of a classical charge density wave (CDW) pumed 
by both the commensurate and the random potentialaE3 
Using the equation ( fl6"| ) , one obtains the following expres- 
sion for the density of states p: 



p{u) 



K 



^ 2 +nK(m 2 +i(u n )) J . 

(77) 

Various correlation functions can also be computed. In 
particular the on-site (CDW) and bond (BOW) charge 
density. For spinless fermions they read 



Xcdw 
Xbow 



{(cU)(c] Cj )) (78) 
(( C ] +lCi +h.c.)(ct +lCj -+h.c.)) (79) 



In the boson representation the 2fcp part of these cor- 
relation function are related to the cos(20) and sin(20) 



correlation functions 



Xcdw oc (-l) a; (cos(2 ( /)( a ;)) cos(20(O))) = K {] (x)(80) 
Xbow cx (-l) a5 (sin(20(a;))sin(2^(O))> = K ± {x)(81) 

Let us now compute these various quantities using 
the result s of the variational method presented in sec- 
tion III B for each of the three phases. 

a Mott Insulator It corresponds to the Replica Sym- 
metric (RS) phase obtained at weak disorder (d/lo < 
2e -1 / 4 ). Because of the non-zero m the whole MI phase 



is thus incompressible (see equation (46)). The MI phase 
is thus the direct continuation of the non-disordered com- 
mensurate one. In the MI phase, the disorder is too weak 
to be able to overcome the gap. 

In the replica symmetric case, p(u>) can be exp ressed 
in terms of the function / defined by equation ([58|) in the 
form: 



p(uj) 



K. 



K 



f(—ix) — x 2 J 2A 

Where x = uj/lu* and cu* = v/£. To perform the analytic 
continuation in Eq. (]8^), we transform ([58]) into a cubic 
equation for / with coefficients depending on x 2 and A. 
Although this transformation adds two spurious solutions 
that do not satisfy /(0) = 0, it proves extremely useful as 
performing the analytic continuation amounts to solving 



the cubic equation for / with x 



Eq. (|8J) implies 



that the phason density of states in non-zero only when 
/ has a non-zero imaginary part. For A < 2, f(—ix) is 
real for : 



l + A-3 



2/3 



(83) 
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FIG. 10: The behavior of the density of states as a function of 
uj for A = 1/2, 1, 2. For A = 1, 1/2 there is a gap in the density 
of states for ui < lj c (X). ui c decreases with increasing A. For 
A = 2, the gap disappears and the density of states behaves 
as p(lo) ~ uj 2 , i.e. there is a pseudo-gap. The pseudogap 
persists in the whole RSB phase. Note that the maximum in 
the density of states decreases as uj increases indicating the 
transfer of spectral weight to low frequencies. 



As a consequence, in the MI phase, there is a gap in the 
phason density of states for uj < uj c — uj*x c 



p{u) = 



UJ < UJr 



(84) 



The physical interpretation of such form for the density 
of states is obvious: no states are available below the gap. 
Thus, in the Gaussian Variational framework, there are 
no discrete two particle states (i. e. excitons) below the 
gap. For uj — > uj c + we obtain p(oj) ~ (uj — wj 1 / 2 . At 
high frequencies, we obtain p(uS) — > K/2 i.e. the density 
of state goes to a constant. A plot of p(u>) is shown on 
Fig.0. 

In the replica framework, we have the follo wing gen- 
eral e xpressions for K\\ and K± (see equations ( A12 ) and 



K\\ = e - ri2( ^ 0) cosh (h2G{x)^ 
K± = e - fi2<5 (o) sinh (h2G{x) 



(85) 
(86) 



In the replica symmetric case, for fi — ► 0, K = K/h 
fixed one finds: 



h6 ^ = m( 1+lj r ,e T 

The resulting correlation functions are: 



(87) 



K±(x) 



sinh 



i^fi + M 



(89) 



For x — ► oo, one has: K\\{x) — > e 16! ° . This implies 
that (cos 20) = d/£. Therefore, in the Mott insulator, 
Charge Density Wave order is still present, but the order 
is reduced with respect to the pure system in which one 
would have (cos 20) = 1. Even for x — > 0, the CDW 
order of the pure system is not recovered. Another inter- 
esting property is that some BOW order is also induced 
at short distances, although BOW order is not present at 
long distance. The presence of BOW order is due to the 
random phase in 4> induced by disorder. When the ran- 
dom phase is of order tt/2, this implies local BOW order. 
However, the positions at which BOW order is obtained 
are not correlated with each other in the system. This 
explains the exponential decay of BOW order. 

b Anderson glass and Mott glass For d/lo > 1-86 the 
system is in the AG phase. In this phase using ( [74]) and 
the expression for G c one finds that the compressibility 

is identical to the one of the pure system \s = • Such 
result is due to the fact that the Gaussian Variational 
Method does not take into account the renormalization 
of K by disorder. Nevertheless, the Replica Variational 
Approximation gives correctly a non-zero compressibility 
for an Anderson glass. We stress that these resiilts are 
valid independently of the presence and absenceEHI of the 
commensurate potential. 

In the intermediate MG phase, with both RSB and a 
gap that is obtained for 2e~ 1 / 4 < d/lo < 1.86, m ^ 0, 
we obtain a zero compressibility. One would therefore be 
tempted to associate this phase with a Mott Insulator. 
Howeve r, the fo rthcoming calculation of the conductivity 
in Sec. [V A 2 will show that this intermediate phase is 
not a Mott Insulator. 

In the replica symmetry breaking case, the formulas 
@ and @ remain valid. However, the function / that 
must be used in eq. (|8|) corresponds to A = 2 in equation 
( p8| ) . This means that as long as there is a RSB solution 
of the variational equations, there is a pseudogap in the 
phason density of states. The behavior of the density of 
states as a function of u is represented on figure E3. 

In the case with broken replica symmetry, one has: 

K\\ (x) = e - 2ri < G (°)> cosh (2h(G(x))) (90) 
K x {x) = e- 2hi > G (°» sinh (2h(G(x))) (91) 

Where we have taken into account the fact that as Ti — > 0, 
TiG c (x) — > 0, and (G(x)) — f Q duG{x,u). Using the one 
step expressions, we obtain: 



17 



H(G(x)) = 



2/ 



e 2, o 



1 - e 1 ^ 



2(1 -M 2 M) 



(92) 



We see that in the Mott Glass phase the CDW or- 
der is still present in analogy with the Mott Insula- 
tor. Such behavior is in agreement with the predictions 
from atomic limit of Sec. Ill A. This time, (cos2^>) = 
cxp[— e*/ (2/i 3 (</>))]. When the system becomes an Ander- 
son Insulator, (cos 2(f)) = which seems to indicate a first 
order transition. Such first order transition is likely to be 
only an artifact of the variational approach. Some sub- 
dominant BOW correlations are also present in the sys- 
tem. They decay exponentially with x and since fi(<p) < 1 
the correlation length of BOW and CDW fluctuations 
is 2lo/n(<p). It is interesting to note that at the Mott 
Insulator-Mott Glass transition, the correlation length is 
continuous. However, there could be a slope discontinu- 
ity which would be characteristic of a second order phase 
transition. 



2 Transport properties 

To differentiate between a Mott and an Anderson glass, 
a crucial physical quantity is the ac conductivity. In the 
Mott insulator, the ac conductivity is zero for frequen- 
cies smaller than the gap whereas in the Anderson glass 
the ac conductivity behaves as u(ui) = w 2 (lnw) 2 in one 



dimensionOtS. Within the GVM, in order to compute 
the conductivity it is sufficient to know m, Si and the an- 
alytical continuation of I(ui n ) to real frequencies- Using 
the Kubo formula, it is straightforward to showcEl that: 



cr(w) 



vK 



7r irK(m 2 + I(— ioj)) — uj 2 



(93) 



where I(iui) represents the analytic continuation of I{u) n ) 
to real frequencies. Introducing the function /, defined 
in (58) one has: 



<j{uj) 



vK 



TTLO* (1 + f{—ix) — X 2 ) 



(94) 



where x = lu/lu*. Similarly to the density of states, the 
behavior of the conductivity is therefore controlled by 
A = |(f ) 3 - One can explicitly check that ( |93| ) satisfies 
the sum rule 



f cLj<j(uj) = 

Jo k 



(95) 



a Mott insulator Let us begin with the conductivity 
for A < 2 i.e. in the Mott Insulator. It is easily seen that 
in order to obtain a non-zero real part of the conductivity 
one must have $sf(ix) ^ 0. As a consequence, the real 




FIG. 11: The real part of the frequency dependent conduc- 
tivity in the MI for A = 1 as a function of frequency. 



part of the frequency dependent conductivity is zero for 
lu < ui c where ui c is the threshold below which the two 
particle density of states is zero (see (|84|). Physically, this 
means that there are no available two particle excitations 
to absorb energy if u> < uj c i.e. at energies below the Mott 
gap. For x > x c (uj > lo c ), the analytic continuation of 
/ to imaginary x has a non-zero imaginary part that 
leads to a non-zero real part of the frequency dependent 
conductivity. For x close to the threshold, 



2 /A 

7! U 



1/3 



Vx^ 



(96) 



As a consequence, for uj > lu c and close to the threshold, 
the real part of the conductivity behaves as dia(uj) ~ (uj— 
loc) 1 / 2 i.e. it is controlled by the available two-particle 
density of states. At large frequency, it can be shown that 
$io~(uj) ~ This behavior can be recovered by a simple 
perturbative calculation in disorder strength. Obviously, 
the conductivity shows a maximum at a frequency ui m — 
(jj*x m (\). The typical behavior of the real part of the 
conductivity for A = 1 is represented on Figure 0. The 
behavior of the threshold frequency ui c a function of A is 
represented on Figure [l^. 

b Mott glass and Anderson glass For A — > 2 the gap 
goes to zero as 2 — A. Quite remarkably, for A = 2, there 
is no gap in the real part of the conductivity although 
the system is still incompressible. The real part of the 



conductivity goes to zero as u) — > as 3t<j(u>) 



. The 



behavior of the conductivity for A = 2 is represented on 
Figure [l3|. As for A < 2 when x — > 00, the real part of 
the conductivity decreases as In fact, this form of the 
conductivity is the one that is obtained in the Anderson 
glass, phase in the absence of any commensurate poten- 
tial^. Moreover, in the GVM framework, the Anderson 
glass is a RSB phaseo. It can be easily seen that for all 




FIG. 12: The variation of the gap in the frequency dependent 
conductivity as a function of A. The gap in the conductivity 
goes to zero linearly for A — > 2. 




(0 



FIG. 13: The real part of the frequency dependent conduc- 
tivity in the MG and AG for A = 2 as a function of frequency. 
For small u, ff(w) ~ uj 2 - 

A > 2, i.e. in all the RSB phases, the scaled conductivity 
is equal to the one of the Anderson glass. The conductiv- 
ity in the MG and AG phases is thus also the one shown 
on Fig. [l3| This remarkable pinning of the scaled con- 
ductivity at A = 2 is a consequence of the marginality 
condition. 



3 General Phase diagram in d = 1 

We have thus identified generically three phases for a 
disordered commensurate system. The bosonization rep- 
resentation being quite general in d = 1 this also applied 
to bosons or spin chains. 

All the previous results having been obtained in the 
limit where K is small, an important question is the range 
of stability of these three phases. Although, in principle 
the variational method could help answering this ques- 
tions, we do not attempt this complicated calculation 
here, and instead give physical arguments. 

It is clear that repulsive enough and finite range in- 
teractions are needed for the existence of the MG. A 
general argument is given in the following section. We 
note here that the case of infinite range (Coulomb) in- 
teraction is a (rather peculiar) example of MG. Indeed 
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FIG. 14: Phase diagram of a one dimensional system with 
both forward and backward scattering random potential. The 
dashed lines correspond to phase boundaries between the 
Mott glass (MG) and the Mott Insulator (MI), the Ander- 
son Insulator (AI) and the Luttinger liquid (LL) phase. The 
separation between the MG and the MI phase in the presence 
of forward scattering disorder is drawn with question marks 
since we do not know how forward scattering affects the com- 
petition of MI and MG phases. 

the one dimensional Wigner Crystal] has a compress- 

2 

ibility Xs = lim^o qi i Q g(i/ g ) = 0, nevertheless it has 
only a pseudogap in the conductivityMS a(ui) ~ u) a . 
One can also show that in a non-interacting system, the 
compressibility gap is equal to the gap for single particle 
excitations. In particular, this means that the interme- 
diate phase cannot exist for K = 1. This result is in 
agreement with the Self Consistent-Born Approximation 
calculation of Mori and Fukuyamafj for the non interact- 
ing case , which do not show any intermediate phase. 
Thus it can only exist for K < K c < 1. 

Let us now give a schematic phase phase diagram, 
which summarizes the effects of both the backward and 
forward scattering in one dimension. As shown in sec- 
tion [A] forward scattering can also lead to gap closure. 
The phase diagram, as function of the LL parameter K 
and the strength of the forward Df and backward D\, 
scattering is represented in Fig. |lj 

B. General arguments and higher dimensions d > 1 

1 Interacting fermionic systems: excitonic argument 

The physics leading to the MG phase is quite general 
and persists in higher dimension as well, as can be under- 
stood through a physical argument. Let us consider the 
atomic limit, where the hopping is zero. One can com- 
pute in this limit the gaps to create both single particle 
and particle-hole excitations (see fig. |l^) Let us consider 
for example fermions with spin with both an onsite re- 
pulsion U and a nearest neighbor repulsion V, with one 
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FIG. 15: Possible excitations in the atomic limit. Drawings 
are made for a chain for clarity but the arguments are valid 
in arbitrary dimension d. (a) Energy cost to add one particle, 
(b) Generic Particle hole excitation, (c) Exciton, where the 
particle and the hole are on neighboring sites. In the presence 
of disorder the gap for excitonic excitations will close first 
leading to the absence of a gap in the optical conductivity, 
but still to an incompressible system (see text), leading to 
the Mott Glass. 

particle per site. Such a system is described by 

H = U fiiffiji + V fijUj 
1 (i,j) 

+ ^W i n i -/x^n i (97) 

i i 

where m — + , and Wi is the disorder potential. 
The energy to add E + i or remove E-\ a particle at/from 
site i are 

E +lii = E^+U + zV-fi + Wi (98) 
E- X ,i = E gs -zV + f i-W l (99) 

and E gs the energy of the ground state of the system 
with one particle per site. If one considers a particle hole 
excitation where the particle moves from the site i to site 
j, the energy cost is S+ij — if i,j are not nearest 

neighbors. On the other hand if the particles are nearest 
neighbors (excitonic excitation), this costs 

\h,ij = U — V + Wj — Wi (100) 

For the pure case, one thus sees from (|8|) and (|l00| ) that 
the gap for creating a single particle excitation is larger 
than for particle hole 

A P = | (101) 
Acx = U-V (102) 

this is the well known excitonic binding that occurs in 
systems with a gap. 

In the presence of disorder one can minimize the sin- 
gle particle gap by choosing the site where the disorder 
potential is minimum, giving 



where we choose Wi — for convenience. On the other 
hand the minimal particle hole interaction corresponds 
to choosing the nearest neighbor pair for which the 
difference in disorder potential is minimal 

A ex = U -V -min\Wj -Wi\ (104) 

(id) 

For an uncorrelated bounded disorder one has 

min(Wi) ~ -W (105) 
min(Wj - Wi) ~ -2W (106) 

Thus, in presence of a nearest neighbor interaction V, 
the particle hole gap closes faster, at W c = (U — V)/2, 
when disorder increases, than the single particle one. For 
an homogeneous system this would simply signals an in- 
stability of the ground state. For the disordered one this 
need not be so, since only a fraction of the sites have their 
gap closing. Thus, in the presence of a small kinetic en- 
ergy the conductivity gap would close near this point, the 
compressibility remaining zero. Within this zero kinetic 
energy model one thus already finds three phases. The 
phase for which the particle hole gap has closed for some 
sites but the single particle gap is still finite can of course 
be identified with the Mott glass. 

Thus the physics of the Mott glass, that has been de- 
rived for finite kinetic energy by the methods of the pre- 
vious sections has its origin in excitonic effects. This is 
quite general and does not rely on any special one di- 
mensional features. One dimension was thus here only a 
tool allowing us to perform the calculation. We thus ex- 
pect the Mott glass to be present in arbitrary dimension, 
and it would be interesting to check either through nu- 
merical calculations or mean field methods whether one 
can recover the properties that we have identified here. 
The excitonic argument also shows clearly that some fi- 
nite range interaction is needed for the MG to appear. 
For a simple Hubbard model both the single particle and 
particle-hole gap would close simultaneously (up to the 
distribution of disorder) and most likely the MG phase 
does not exist. In the presence of finite range interactions 
the MG glass can be stabilized. 

A similar construction can be made for the spinless 
case, although it involves longer range (third neighbor) 
interactions. 

According to this physical picture of the MG, the low 
frequency behavior of conductivity is dominated by ex- 
citons (involving neighboring sites). This is at variance 
from the AG where the particle and the hole are created 
on distant sites. This has consequences on the precise 
low frequency form of the conductivity such as logarith- 
mic corrections. In addition since the excitons are neutral 
objects, although they can participate to the optical ab- 
sorption, they need to be broken to give a d.c. current. 
One can thus give a naive estimate of the conductivity in 
the MG 

a^n ex e~^ (107) 
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FIG. 16: Phase diagram of bosons in the limit T — > 0, 
h -> 0, W is the strength of disorder and g the amplitude 
of the commensurate potential for a fixed repulsive interac- 
tions. For weak disorder the boson system should retain per- 
fect topological order (i.e. no defects such as dislocations). 
Whether the line for which topological order is lost (dashed 
line) enters the MG phase or not is an open question. It is 
represented here for d = 3. The incompressible phases are the 
Mott Insulator (MI) or Mott Glass (MG). The compressible 
ones are either the Bragg Bose Glass with perfect topological 
order (BBG), or the Bose Glass (BG) if topological order is 
lost. 



where n ex is the number of excitons in the ground state 
and V is the typical excitonic binding energy, which de- 
pends only weakly on the disorder. 



2 Consequences for other systems 

The above arguments also directly apply to other sys- 
tems. In one dimension the spinless fermions can be 
mapped to a disordered spin systems. In that case the 
commensurate phase can either come from an antifcr- 
romagnetic staggered field, or more reasonably from a 
spin-Peierls distortion of the lattice. Such a perturba- 
tion would force the spin to lock into a singlet state. The 
disorder would be a random magnetic field. 

Another system of interest is provided by hard core 
bosons. In one dimension, one can use exaGtlypthe phase 
Hamiltonian to represent interacti ng boso nalja, but the 



excitonic arguments given in Sec. IV Bl would also ap- 
ply to interacting bosons in higher dimensions as Well Ty 
that case the Anderson glass becomes the Bose glassEBo. 
For classical systems the phase diagram is shown on 
Fig. [h]. Even in the Bose glass phase, for weak disor- 
der perfect topological order (stability to dislocations in 
the lattice) r caa. persist in d = 3 resulting in a Bragg Bose 
glass phaseoEl On the other hand in any dimension the 
Mott insulator should exhibit perfect topological order. 
Thus an interesting and open issue is whether this topo- 
logical order also subsist in at least a portion of the Mott 
glass phase. 



V. CONCLUSION 

In this paper, we have investigated the competition 
of a random and a commensurate potential. This ques- 
tion is relevant for host of physical systems, ranging from 
one dimensional interacting fermions or bosons, to clas- 
sical systems in the presence of correlated disorder. The 
commensurate potential induces an incompressible Mott 
insulating phase with a gap in the conductivity. On the 
other hand, disorder induces a compressible Anderson in- 
sulating phase. While naive expectations predict a direct 
transition between these two phases, we find that if in- 
teractions are repulsive enough, an intermediate phase, 
the Mott glass, does exist. Although this phase is incom- 
pressible, as a Mott insulator, it does not have a gap in 
the optical conductivity, in a way similar to the Anderson 
insulator. 

To obtain this phase we had to go beyond standard 
renormalization group techniques which are perturbative 
in both the commensurate and disorder potential. We 
therefore have used bosonization associated with several 
non perturbative techniques. The first one is a replica 
variational method, that allows for a complete calculation 
of the various physical observables such as the conductiv- 
ity. The second method is a functional renormalization 
group, which is perturbative in d = 4 — e dimensions and 
is well suited to study the transition from the Mott glass 
to the Mott insulator, as well as equivalent classical sys- 
tems. In addition we have looked at the limit of zero 
kinetic energy, both for the bosonized Hamiltonian, and 
directly on the fermion problem (both for the spinless 
problem and for the problem with spin) . The later yields 
a very general argument in favor of the existence of the 
Mott glass in any dimension. It also shows that the un- 
derlying mechanism for this phase is the creation of low 
energy bound states (excitons) coming from the competi- 
tion between interactions and disorder. These excitations 
play no role in the compressibility but contribute to the 
optical conductivity. 

This phase could be observable in systems close to a 
metal insulator transition, such as oxides, provided that 
one can measure simultaneously the optical conductiv- 
ity and the compressibility. Numerical simulations on 
disordered boson systems could be prime candidates to 
observe this effect. Note that since all the phase have a 
finite correlation length, this should be observable even in 
moderately small systems. Many problems remain open. 
In particular, it would be interesting to understand in de- 
tail the effect of a chemical potential on the Mott-Glass 
phase. Another open problem is the effect of temperature 
on the Mott Glass phase. Finally, it would be interest- 
ing to investigate the possibility of aging dynamics in the 
Mott Glass. 
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scattering. For fermions, such an approximation is surely 
justified when 3/2 < K < 2. Then, the backward compo- 
nent of disorder is irrelevant and can be neglected. In the 
other case, K < 3/2, backward scattering will be relevant 
and drive the system into an Anderson glass state. 



1. solution of the variational equations 

The action of the problem is: 



APPENDIX A: FORWARD SCATTERING 
DISORDER AND PERIODIC POTENTIAL 

In this appendix we examine the effects of the forward 
scattering disorder r\ and neglect altogether backward 
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We use the GVM ansatz (j40D with 

vG-£(q,u>) = (M ^\ fc - a ab (q,uj) (A3) 

as in the case of the backward scattering disorder. Using 
© and (||) the variational energy F var for the forward 
scattering problem is: 
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Minimizing (A4) with respect to G(q,u>) gives the varia- 
tional equations 
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It is easy to check that the equations ( |A5| ) only have 
replica symmetric solutions in contrast to the case of 
backward scattering. 

Using the standard techniques for inversion of matrices 
in the limit n — > 0c3, one finds the following expressions 

for G c = G a a + X/b^a G a b- 



G c (q,u n ) 



J^(u* + {vq) 2 )+m 2 
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and for G(q,u) n ,u): 
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FIG. 17: The reduction of the 
disorder strength 
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where: 
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and £ 2 = v 2 /(irKm 2 ). One has: 
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Leading to the self-consistent equation for £ 2 
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Where we have defined ^ 1 
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It is straightforward to show that (A10) has two solu- 
tions for h/d > e/2 and no solutions otherwise. In the 
first case, the solution with < 1/2 is a spurious solu- 
tion, as can be seen by taking the limit of zero disorder 
(h — > oo). Physically, the Mott gap is preserved (but 
reduced) as long as li/d > e/2, whereas for l\/d < 1/2, 
the Mott gap disappears. The transition gapped-gapless 
appears first order in the GVM, which is likely to be an 
artifact of the variational method. Note that the condi- 
tion for the transition 1%/d ~ e/2 is in agreement with 
strong coupling extrapolations of the perturbative RG 
treatment. 



2. Correlation functions 

Another advantage of the GVM is to allow for the cal- 
culation of the correlation functions. Using ([73]) one sees 
that the density-density correlation functions is given by 

/J' i . e i20(a:,T) e ±j20(O,O)\ = g-fi[2G(0,0)±2G(x,"r)] (All) 

Which leads to: 



K x (x,t) = (T r sin(2(/)(x,r))sin(2^(0,0))) 

= e- ri2(5 (°<°)smh?i2G(:z;,T) (A12) 

Aj|(x,r) = Wt cos(2<f>{x, r)) cos(20(O, 0))) 

= e- h2€;(0 ^coshh2G{x,T) (A13) 
Since G(x, t) = G c {x, t) + G(x) and lim^o hG c (x, r) 
0. only the static correlations survive. It is straightfor- 
ward to show that: 
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Leading to the following expressions for the correlation 
functions: 



K±(x) = e 1 sinh 



Kii(x) = e *I cosh 
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(A15) 
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-.1 

It is easily seen that for £ < oo, lim^^oo K» — e 1 . In 
the absence of disorder, this limit would be exactly one. 
Forward scattering disorder thus leads to a reduction of 
the Charge Density Wave long range order. For £ — * oo, 

one recovers K m (x) ~ e 'i , a result that could have been 
derived directly. 



APPENDIX B: SADDLE POINT EQUATIONS 

We derive in this appendix the saddle point equations 
obtained by minimizing the variational free energy (Ej). 
Using (ph and <M) we get 
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Varying in (Bl) with respect to G we get the following 
saddle point equations: 
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with G c = G aa + Ysb^ta G °-b and 

B a b(x, t) = G aa (x, t) - G ab {x, t) 



(B4) 



As was the case for fermions in a random potentia 
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expression of the replica off diagonal self energy: 
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leading to the following simplified (B2) and 



We still need to perform the analytical cont inuation from 
positive integer n to n = in equations (B2) and (B3). In 
the GVM this is done assuming that for n — > 0, the G a b 
become hierarchical matrices. Using the Parisi paramcL- 
tcrization of hierarchical matrices in the n — > limitca. 
give the following equations: 
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where u £ [0, 1] is the Parisi parameter replacing the 
discrete replica index a. 



APPENDIX C: SOLUTION OF RSB 
EQUATIONS 

We want to solve the RSB saddle point equations 
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Using the inversion formulas for hierarchical matricesE 2 ], 
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Here, we have assumed that when /? goes to infinity, m c 
goes to zero in such a way that /3u c = 5 remains finite. 
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Next, we derive a self-consistent equation for m by taking 
the h — ► limit of the equation: 
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We use first the general inversion formulae 2 ]: 
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In which we have: 
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leading to the expression for m in the (3 — > oo limit: 
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A final equation for the breakpoint u c is needed to close 
the system of equations. As was discussed in Ref. ||(] the 
physical choice corresponds to the so called marginality 
of the replicon condition which yields to I(ui n ) oc| | 
and to : 



(7tq) 2 (to 2 + S1) 3 / 2 
Using the quantities 
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Where l and d are defined respectively by Eqs. (J51 
and (|52]). The reduced variable \x is defined in such way 
that The point at which the replica symmetric solution 
becomes unstable has [/, = 1. 

the self-consistent equations are rewritten: 
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To solve (C24) we introduce ip = r] 1 ^-^. Physical solu- 
tions have 77 > 0, < \x < 1 an d thus ip < 0. Excluding 
the solution /i = 1 from (C24) we obtain the following 
equations in terms of <p: 
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We have thus reduced the self-consistent equatio ns to a 
s ingle equation for (p. Putting <p — in Eqs. (|C2 7|)- 
( |C28| ) we obtain /i = 1 and we recover the condition (|60|) 
on d/lo, i.e. the limit of validity of the RS solution. A 
plot of F(ip) in shown in Fig. As can be seen from 
this plot, F has a minimum for ip = <p c . T his i mplies 



that there can be no solution of Eqs. ( C27 )-( C28 ) when 
lo/d < \J F(<p c )/2. The physical values of <p are thus lo- 
cated in the interval [<p c , 0]. Numerically, it is found that 
<p c = -3.4325 ±0.0025 and F(tp c ) = 1.15338. The corre- 
sponding critical value of lo/d is then lo/d = 0.536977 i. 
e. d/l = 1.86.... 



km) 




FIG. 18: The graph of F(cp). Only the region with cp < 
is physical. F has a minimum for ip = (p c . When 4(lo/d) 2 is 
smaller than F(tp c ), equation (C28) has no solution. 



APPENDIX D: FUNCTIONAL 
RENORMALIZATION GROUP APPROACH 

In this Appendix we detail the analysis using the Func- 
tional Rcnormalization Group method of the effective 
model ( |66|) in presence of a mass term and correlated 
disorder in dimension d. We use the notations of the 
classical equivalent model (|22]). The method is a Wil- 
son momentum shell integration which it is an extension 
of Ref. [32] to the case of a finite mass m 2 > 0. Similar 
extensions can also be found in Ref. |5^ , p5| , p6U6l] . 

We start by studying T r j = 0. In the quantum problem 
this corresponds to the limit K — * 0, Ti — * 0, K = K/% 
fixed (see Eq. (^)). We consider the ground state which 
is t independent (p{x, r) = (p{x). It is in this limit that 
the GVM method revealed the presence of the Mott Glass 
phase. It is more convenient to work with the function 
A(0) = —R"((p), where the bare R(<p) has been defined 
through (^|) . One first defines the running dimensionless 
disorder: 



A,(0) = ^Af-%(0 



(Dl) 



with A d — Sd/(2-K) d , which is found to obey the T ct = 
FRG equation: 



diA = eA-f(l) (A" 2 
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A - A(0) (D2) 
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The UV cutoff is reduced to A; = Ae 1 (A ~ 1 /a where a 
is lattice constant). One can check that for m = 0, (D2) 
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reduce to the one derived in Ref. ^2|. This equation turns 
out to be identical to the one describing point disorder 
in dimension d. .— . 

It is well knowr£3 that in the case m = 0atT = 0a 
cusp develops at the origin for I — + oo. One finds that 
|A(0,oo) - A(0,oo)| oc \<f>\ for <f> -> 0. This implies that 
lim/^ +oc A" (0,1) = — oo. It is thus important as a first 
step to analyze the cusp generation in the c ase m ^ 0. If 
we define A 2 (l) = -A"(0, 1), we have from (p2|): 



d t A 2 = eA 2 + f(l)A 



This differential equation has for solution: 



dx 



A 2 (0) A 2 (l) 2 
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Where A2(0) is the bare disorder. Introducing the Larkin 
length R c in the absence of a mass (fi — 0), defined as 
the length scale at which A 2 diverges i.e. 



A 2 (0) 



1 [ {R c /a) dx 
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One also obtains an equation that determines the Larkin 
length in the presence of a mass R c (n) defined as the 
length where A 2 = +oo for a non zero /j as a function of 
R c and fi : 



(flc/a) 2 



dx 




4-d N 
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(D7) 

(D7) has two types of solutions, one with R c = oo for 
weak disorder and another one with R c < oo for stronger 
disorder . As discussed in the text, this means that there 
are two phases, one in which disorder is strong enough to 
generate a cusp and a second one in which the flow is cut 
by the presence of the mass before a cusp can be gener- 
ated. The former corresponds to the Mott Glass phase, 
while the second one corresponds to the Mott insulator 
phase. The equation of the transition line betw een these 
two phases is obtained by setting R c = oo in (D7) and 
reads R c = R* c (fi). At small /i and for d < 4, we find 
that it behaves as: 



R*c(p) c(d) 



(D8) 



where C(d) is a dimension dependent constant. 

The physical quantity which is directly affected by the 
presence of the cusp is the tilt modulus ca(l). One finds 
that it satisfies the RG flow equation: 



ftlnc44(0 = -A a (0/(0 



(D9) 



while c remains unrenormalized. So clearly, either A 2 (l) 
diverges sufficiently fast and c^(l = +oo) = +oo (in- 
terpreted as the Mott glass) or the mass cuts off the di- 
vergence early enough and c^(l = +00) remains finite 



(Mott insulator). ^,From the FRG it is possible to com- 
pute exactly the large / behavior of the tilt modulus. For 
that, and to make further progress in the analysis of the 
two phases, we need to first consider the full flow of the 
function Ai(<f>). It can be shown easily that the solu- 
tion of the flow equations at fi ^ can be obtained as a 
function of the solution at fi = in the following way: 
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,l) = h(l)A ll = (<t>,t(l)) 
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l + etidl> Tl 
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with the same initial condition. The behavior at large I 
is the following: 



h(l) 
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where R*((J<) — ae'«W W as defined above. From ( |D9| ) 
it is then easy to see that 044(1) = 044(0) A2^)/A2(0)e d 
which yields e.g. C44(+oo) in the no cusp phase as: 



c 44 (+oo) 



'Rc 
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C44 (0) (Rcy fMtiy 
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One thus finds that the renormalized tilt modulus di- 
verges as one approaches the transition as: 



c 44 (+oo) - (R c - R*(V)Y 



In d = 4 one has instead: 



c 44 (+oo) = c 44 (0) 
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In all cases one has 044(00) — » +00 in the cusp phase. 
On the contrary in the no cusp phase 044(00) remains 
finite. We expect that having 044(00) — > +00 leads to 
no conductivity gap but having 044(00) < 00 produces a 
conductivity gap. In the cusp pJaase one can also expect 
that a term \d z u\ is generated^. Such term give rise 
to the transverse Meissner Effect. The critical field h c i 
needed to bend vortices can be easily computed from the 
FRG. 

We are now in position to estimate correlation func- 
tions in the case of point disorder or their z (i.e r) in- 
dependent part in the case of correlated disorder (i.e at 
co n = 0). One has: 



M?M-?)> = r(?) - e M T(qe l , A(l), me 21 ) (D16) 

In the regime qa ~ 1, the correlation function T can 
be obtained by perturbation theory in A. Our strategy 
to obtain correlation functions^ is therefore to integrate 
the RG equations until qae 1 ~ 1. At this point, we can 
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calculate the correlation function T perturbatively and 
deduce f. We obtain: 



f(?) 



A(0,Z = ln(l/og)) 



{aq) d ^ i {q 2 + n 
So that for d — 2 one gets: 
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Note that the static two point correlation functions or 
equivalently the correlations for point disorder do not 
exhibit a sharp transition. 

It is crucial to check that the transition we found for 
T c i = (i.e h = 0) survives at hnite temperature (finite 
h). In the original quantum problem this corresponds 
to K > i.e. whether the intermediate phase exists for 
interactions that are not infinitely repulsive. The RG 
can be performed at hnite T c /. Introducing the effective 
running temperature: 



f x =Td-^±=k(l) 
2^1(0 

fc(Z) = (l + M e 2 'r 1/2 
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one finds that the FRG equation becomes: 

d,A = eA + T,A" - /(/) [(A') 2 + A" [A - A(QMl) 



Note that in the quantum parameters To ~ K. ._. 

In the absence of a mass, fj, = 0, it is easy to seeEj that 
the temperature 7} runs to exactly zero at a finite length 
scale I* (f Q ) - l c ~ c(d)f for small f with l c = ln(R c /a) 
the Larkuruacale. This is because for I > l c the cusp is 
roundedE§N at finite Tj with : 



. m A*'(0+) 2 2 1 
J- 1 i 



(D22) 



where A*(</>) is the T = fixed point function, and \ a 
numerical constant. Thus one can write: 

ftlnf, = l-d-^lnc z (0 (D23) 



T, 



(D24) 



This yield that d{Ti « — xe 2 and thus the temperature 
vanishes beyond the scale I* (To) — l c ~ -^jTq. 

It is thus clear that if /i is small enough so that 
I* [To) "C introduced above, the temperature will 

vanish before the term f(l) starts deviating from 1 and 
change the behavior of the solutions. Thus at small non 
zero temperature the divergence of C44 is not suppressed 
and the transition survives. 



A more detailed analytical study can be performed 
noticing that the relation between the solution at finite 
[i and zero mass: 



(0,O = /*(OA U=O . # ,(<M(O) 



(D25) 



with the same functions h(l) and t[l) as above and T; = 
Ti/(h(l)f(l)). It confirms the above conclusions but will 
not be detailed here. 

Note finally that the above RG procedure uses that the 
thickness L is constant. Since h runs to zero it means that 
j3 runs to infinity, and thus that the temperature is also 
irrelevant in the quantum system. 
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